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lo  Introduction 


Previous  technical  reports  have  dealt  wi  th  light  scattering  by  a single 
spherical  particle,  or  by  a raono-disporoe  system  of  such  particles.  T/hen  tie 
scattering  by  a hetero-di sparse  system  of  parti  cles  is  used  to  determine  the 
particle  s1k©  distribution,  vre  are  ooncomcd  with  the  determination  of  an  unknown 
function.  To  detomine  such  a function  completely,  it  is  evidently  neoeBsary  to 
observe  sone  quantity  which  depends  on  a parameter  which  can  be  varied  continuously 
Two  such  parameters  imediately  suggest  themselves i the  wave-length  of  the  light 
and  the  angle  of  observation.  In  this  report  only  the  first  of  these  is  considered 
and  attention  is  focussed  on  the  question  of  what  can.be  determined  as  to  particle 
site  distribution  from  turbidity  measurements  whon  the  wave-length  is  varied. 

In  is  assumed  throughout  tl»t  the  scattering  particles  aro  all  spherical  in  shape 
and  of  tho  sane  refractive  index  (independent  of  the  wave-length),  and  that  the 
concentration  is  sufficiently  low  for  multiple  scattering  to  -be  neglected, 

2„  Tho  internal  equation  for  the  distribution  function  and  its  solution. 

Tho  scattering  cross-section  for  a aiiiglvi  spherical  particle  of  radius  X 
when  illuminated  by  monochromatic  light  of  wave  -length  \ can  be  written  in  the 
form 


cr  - 


at  t ' ’ 


where  k * 2vA, 

n 

and  f is  a funotion  which  'ms  boon  extensively  tabulated."  It  depends  on  the 


■^Seo,  for  instance.  Technical  Reports  1,  Z,  5 of  this  soriec;  Tables  of 
Scattering  Functions  f or  Spherical  farticleB.  National  Bureau  of  Standards 
^Wasliington,  1943)j  I?.  0.  Gurprecht  and  C,  II.  Gliepoevitch,  Tables  of  lights 
scattering  FuncMonr  for  Sphe rical  Particles,  (U.iivor  sity  of  Tliciiigan,  1'151)0 
The  function  f is  .identical  wi-ch  the  function  2 of  Technical  hoport-  ??c,  1. 

In  terns  of  the  e nattering  coefficient  X,  j ~ ^ [<(*(/) 

whore  d~R.x.^  ' ' s > 


rofractive  iv  v<  as  ~/eli  ac  on  1 c:;  we  h'-'.crv,r,  or't:  tho  depenoonce  on  m. 

let  the  rm‘ibo7‘  of  scattering  particle"  per  writ  volute  with  radii  in  the 
range  ( T,  , t-uix  } be  n(  X ) J- .t.  , and  lot  t!»  turbidity  at  wave-length  \ bo 
denoted  by'f^X),.  Neglecting  multi j a scattering,  vie  then  have 


Fft'-l  = 


(2) 


where 


A 


(3) 


The  problem  is,  tlisn,  to  colvo  the  integral  equation  (2)  for  the  function  "Wpx}* 
the  function  F(!A  and  the  kernel  being  regarded  as  known. 


KX) 


If  we  suppose  that  F(k)  is  known  for  all  values  of  k from  0 to  oO  (i.o« 


is  known  for  all  wave-lengths  from  0 to  CO  ),  then  (2)  is  of  a type  for  will  oh 

O 

a rigorous  solution  in  possible  by  fjollin  transformers • **  Unfortunately,  however, 
this  solution  is  of  little  practical  use  in.  tho  present  caoo.  For  in  practice 
the  turbidity,  and  lser.ee  F(k),  will  only  bo  known  for  quite  a limited  range  of 
a range  not  sufficient  to  deterrti.no  the  TTellir.  transform.  ox’  F(ro)<>  Furthor,  even 
if  F(  it.  ) 7,'orj  known  for  a sufficiently  wj.de  range  of  k,  it  would  then  be 
necessary  to  take  account  of  the  variation  of  refractive  index  with  wave-length, 
and  this  would  destroy  the  special  A’oru  of  the  equation  (2).  A further  objection 
is  that  tho  Hellin  transf or:  i so?:  tin:  • • ild  require  fairly  extensive  numerical 
computation,,  , 

Various  other  net-hods  of  solving  (2)  ware  therefore  cousidorrxi , ,Tn  one  of 
these,  it  war,  assumed  that  tho  distribution  was  narrow,  about  soma  iioan  radius  r'  0 


^Eee,  for  instance.  F»  do  litchaarch.  Introduction  to  the  Theory  ox*  Fourier 
Integrals,  (Oxford,  19*8),  316* 


say,  and  e Taylcr  . oriec  for  tea  fanatic...  .^Jtc'  :'i;  neighborhood  of  )bt0  v/as 
used,  the  derivatives  of  tho  function  f being  found  uunorlcally,  The  calculation 
of  the  higher  derivatives  proved  unse.ti. afechory,  however,  owing  to  irregular 
differences  in  thn  tabulated  function. 

Further,  this  netiiod  would  be  of  no  uso  for  a broad  distribution,  and  was 
therefore  abandoned. 

Another  method  attempted  vm.n  to  replace  tbs  integral  equation  (2)  by  a set 
of  linear  equations,  1110  requisite  matrix  inversion  was  performed  by  the  VJayne 
Computation  Center,  T.'tis  met'ood  also  proved  unsatisfactory,  due  to  the  fact  that 
some  of  the  eloTiients  of  the  inverse  matrix  vrere  very  large  compared  with  others, 
30  tliat,  as  the  funotion  F(k)  in  (2)  is,  in  practice,  subject  to  experimental 
errors,  the  results  deduced  in  this  way  would  be  unreliable , This  behavior  of 
the  inverse  matrix  io  evidently  connected  with  tlie  fact  t!»t  the  function  f is  a 
rapidly  increasing  function  of  its  argument  (see  on).. 

The  method  finally  adopted  was  the  following,  Although  it  has  not  yet  been 
brought  to  a numerical  conclusion,  it  is  bei:  evert  that  it  should  give  reliable 
results  subject  to  certain  restrict!  on.n  on  the  distribution  function.  In  the 
first  place,  we  shall  suppose  that  the  particle  radii  lie  within  a certain  range, 
say 

x0i-A  x,,  (4) 

and  tliat  observations  of  turbidity  are  nr.de  at  wave-lengths  also  lying  within 
a certain  range,  say 

kc-Ah4h4k-i-At_<-  <6> 

We  must  then  write  (2)  in  tho  form- - ■— 

rxc1-&X. 

F(t')=  f(.Mo,U)clx,  <6> 

Jx~£>x. 

& 

where  kQ  — --  4 4 ^ 4 ^ ' 


Let  ua  .'to'"  rz.h 


^Institutions  a 
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xo+  £,x-  y/ 


L+-^k-  b/  • 


then  (6)  can  be  written* 


where 


4 W)  - f V)  N Mix.'' , H 4 k't  0 
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(8) 

(9) 

CIO) 
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4*  C h/y  x/ ) - 5 ( k x)  ? 

N(x ')  =l  'Hlx)  A x (is) 

We  observe  that  i|(se^)ibe*“  is  ti*>  SwJiiber  oC  particles  par  unit  volume  With  radii  between 

X64-  &X-  x/  Xfc+-  £ X(XV  4s/)  • 

Let  us  row  suppose  that  the  kernel  'fCk/JX/)  i«  expan-k-d  in  a series  of 
Legendre  polyror.  iul  s r 

Q& 


(13) 

xi-  ut11  j i-*-  y* 

1 *.,z»0  < <f  y 

Multiplying  (3)  by  P-j  {k  / and  integrating  with  rospeet  to  it1  from  no  10  we  obtain t 


where 


od  * r ! 

ZL  Ci;  N;  aktil  I ^hbP-tfe') 

j-o  J <?  a.  j . : ' ' ■ > 

’#  ^ ' *r*  j 

Nj  ~ j A^  ( *-*)  F (x/)  4 j4. 


(14) 


(is) 


Equations  (14)  now  constitute  -ar.  infinite  not  of  ilnoar  equations  for  the  detons. 

i 

alimtlon  of  the  unknown  cofista >**-5  h-6  living  found  these,  $(x  ) in  given  by 

00  ^ 

M(xx)~  *£_  A>  T*  (*/)  . 

4 ~ * * 


In  practioo.  of  course,,  a fi-dto  ::  dur  of  i:«rmn  In  ( 13  ) „ (14),  (1(5)  must  be  used. 

Any  system  of*  orthogonal  functions  could  be  used  in  place  of  the  Legendre 
polynomials*  Those  particular  functions  are  suggested,  however,  by  tko  foot  that 
the  function  f can  be  well  approximated  over  a fairly  v/ide  range  of  ite  argument 
by  a polynomial  of  fairly  lav/  degree,  say  of  degree  p„  Taking  account  of  (7), 

(8),  (11),  this  means  that  (k\  x^)  can  bo  approximated  by  a polynomial  which 
is  of  degree  p in  each  of  the  variables  k , x , and  hence  by  using  only  Legendre 
polynorriala  of  degree  jC  p in  (13)0  With  this  approximation,  (14)  beocnea  a set 
of  (p  ♦ 1)  equations  in  (p  + 1)  unknowns,  For  the  ranges  of  partiole  size  and 
wave-length  so  far  considered,  p could  bo  quite  a arm  11  number— say  fro®  3 to  5, 
Whether  a rathor  small  number  of  terns  in  the  eerier.  (16)  yiold3  a good  approxi- 
mation to  the  true  distribution  function  depends,  of  course,  on  the  nature  of  tins 
function  and,  specifically,  on  whether  it  cau  be  well  approximated  by  a polynomial 
of  low  degree.  The  trend  of  the  vnluo3  of  the  coefficients  Nj  should  give  a fairly 
good  indication  a3  to  the  rapidity  of  convergence  of  the  series  (16), 

Irrespective  of  the  convergence  of  (16),  however,  a knowledge  of  a certain 
number  of  the  coefficients  Hj  provides  us  with  useful  information  about  Lite  dis- 
tribution function.  For,  from  (7),  (12)  and  (IS),  we  see  that  if  the  first 

(p  + 1)  coefficients  TT.  axe  known,  then  the  first  (p  1)  nonsnto  of  the  distribution 

< J 

with  respect  to  thu  origin,  namely 

• 

N ’ ~ \ ,-jr.  t.  K ) y-  ^ X U - Q \j  ■ ■■  ■ j V>  (17) 

^ J X',"  A X 

can  be  found,  Tiie  moment  v.u  givrs  the  total  number  of  particles,  while  l\j/l lQ  gives 
tho  nean  radius  of  the  particles,  higher  no* .’one a furnish  tho  xasun  square  deviation, 
moan  cubic  deviation,  and  so  on,  so  that  considerable  information  as  to  this  nature 
of  tho  distribution  can  be  obtained  i . fci.i-'  way. 

If  it  is  desired  only  to  calculi- In  tha  me  .ants  Mj,  when  the  function  f is 


approid  mated  by  n p lynonial the  ubent*  procedure  nay  be  ahortonod  as  follows:  let 


G 


us  write,  at;  tti  ia 


4-(  ,1..  c-  (jj  >.y . 

<j  ~°  J 

Then  substituting  this  in  (G)  and  proceeding  ae  befo.ro  s we  obtain 

4 


where 


t 


ci  i (ivr,,k  u'dp  Cv,/V)h'. 

O / ' A,  ^ ✓ 

-'-I 


(IB) 


(19) 


Equations  (18)  are  a sot  of  (y  * 3 } equations  to  determine  tluo  (p  + 1)  moments 
It.,  ?!•. , Mn<  i’he  eoofficis.-n ts  Q_  can  bo  calculated  once  and  for  alio 

Tho  integrals  occurring  in  (it)  (or  can  be  calculated  by  nuner.i cal  integration 


fron  tho  observed  turbidity  or  — poiliaps  rather  noro  simply  * ~ by  approximating 
the  function  •j>t  r;/)by  a polynomial  in  k-^  and  then  expressing  the  various  powers  of 
k1  in  terms  of  I/.'gendre  polynomials.  Since^frora  (19),  a.^  - 0 if  i > j,  the 
equations  (18}  cr-.n  bo  solved  without  difficulty.  An  experimental  check  can  he 
obtained  from  ;ho  value  of  the  moment  Ug.  For  the  total  mass  per  unit  volume  of 
the  scattering  particles  iaij  p Mi,  where  p is  the  density  of  tho  particle  material, 
and  this  quantity  will  be  luiown  or  can  bo  dote  mi  nod. 

In  connection  with  t’lis  method,  tJiO  following  point  slwaiu  be  observed i the 
function  f is  a rapidly  increasing  function  of  its  argument  over  the  useful  range, 
especially  for  small  valuos  of  the  argument  ( J^fir  « ), 

which  means  that  tl*e  larger  p:;i- titles  ooriGributo  much  more  to  the  scattering  than 
the  smaller  onus,  Consequordly  wuiv;  to  the  experimental  errors  in  the  determination 
of  turbidity,  it  v/iil  not  ue  possible  to  do to  nr  ti.no  accurately  tns  distribution 
function  for  tho  oral ler  values  of  x if  tho  distribution  is  a broad  one.  The 
moments  in  should,  uevurtladosc , bo  given  v;i  ;:h  fail  accuracy  r. 

Another  relator!  j o.vr  is  the  following  : although  the  functi  an  f can  he  roll 
approximated  b..  a ^olynonial  of  lev/  degree  over  u useful  rr-.ugo  of  tu  argument. 


7 


3uoh  a rapre39.1tf.tion  gives  a larrp;  rcontage  error  (though  a email  absolute 
error)  for  snn.ll  values  of  the  ar.511r.oftt  (assuming  tint  snail  values  of  the 
argument  have  to  be  used,  which  i s tlie  case  if  the  lower  limit  of  particle  radii 
iB  snail).  This  Yri.ll  not  natter  if  the  distribution  is  a broad  one.  If.  however, 
the  distribution  is  such  tliat  an  appreciable  freotion  of  the  relevant  values  of 
lex  are  snail,  tlje  rosults  obtained  will  be  misleading.  To  take  an  extreme  case,  if 
the  approximation  4"(  U xi)  a const,  /Ji.  x.'f  (Rayleigh  scattering)  is  valid  over  the 
whole  relevant  range  of  valuer  of  kx,  then  evidently  all  that  oan  be  determined 
about  the  distribution  is  the  value  of  the  moment  l!g. 

It  is  lioped  to  issue  a further  report  shortly  where  numerical  data  will  be 
given  which  yjIII  enable  the  method  to  bo  carried  tlirough  explicitly  and  without 
too  much  labor.  It  i o hoped  also  to  present  comparisons  of  particle  size  distribu- 
tions obtained  from  turbidity  measurements  vri.  fch  those  obtained  eleotron-nioroscopi- 
cally. 

3_o  Conclusion. 

It  chould  be  possible  to  deter."..’  no  tho  first  fev/  moments  of  the  particle 
radius  distribution  from  turbidity  ; voasureneiits  wit;',  fair  accuracy',  provided  that 
not  too  wide  limits  can  ha  sat  to  "ho  particle  radii,  ?md  tliat  not  too  large  a 
fraction  of  the  particles  have  ' ’ ' « radii.  From  these  moments  the  distribution 
itself  con  bo  calculated  if  t!>o  aorv  a a of  Legendre  polynomials  (1G)  is  sufficiently 
rapidly  convergent. 


